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Donor states in Si nanodevices can be strongly modified by nearby insulating barriers and metallic 
gates. We report here experimental results indicating a strong reduction in the charging energy of 
isolated As dopants in Si FinFETs, relative to the bulk value. By studying the problem of two 
electrons bound to a shallow donor within the effective mass approach, we find that the measured 
small charging energy may be due to a combined effect of the insulator screening and the proximity 
of metallic gates. 
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I. INTRODUCTION 

For over a decade dopants in Si have constituted 
the key elements in proposals for the implementation 
of a sohd state quantum computer.^"™ Spin or charge 
qubits operate through controlled manipulation (by ap- 
phed electric and magnetic fields) of the donor electron 
bound states. A shallow donor, as P or As in Si, can bind 
one electron in the neutral state, denoted by D°, or two 
electrons in the negatively charged state, denoted by D^ . 
Proposed one and two-qubit gates involve manipulating 
individual electrons or electron pairs bound to donors 
or drawn away towards the interface of Si with a barrier 
materialji?'^ In general, neutral and ionized donor states 
play a role in different stages of the prescribed sequence 
of operations. 

In the proposed quantum computing schemes, donors 
are located very close to interfaces with insulators, sepa- 
rating the Si layer from the control metallic gates. This 
proximity is required in order to perform the manipula- 
tion via electric fields of the donor spin and charge states. 
The presence of boundaries close to donors modifies the 
binding potential experienced by the electrons in a semi- 
conductor. This is a well-known effect in Si MOSFETs,'^'^ 
where the binding energy of electrons is reduced with re- 
spect to the bulk value for distances between the donor 
and the interface smaller than the typical Bohr radius of 
the bound electron wave- function. On the other hand, 
in free-standing Si nanowires with diameters below 10 
nm, the binding energy of donor electrons significantly 
incrcascs'^'^'' leading to a strongly reduced doping effi- 
ciency in the nanowires J^ 

The continuous size reduction of transistors along 
years, with current characteristic channel lengths of tens 
of nanometers, implies that the disorder in the distribu- 
tion of dopants can now determine the performance, in 
particular the transport properties of the devices .J^^-ii In 
specific geometries, like the nonplanar field effect transis- 
tors denoted by FinFETs, ^^ isolated donors can be iden- 



tified and its charge states (neutral 13°, and negatively 
charged D~) studied by transport spectroscopy. 

The existence of D" donor states in semiconductors, 
analogous to the hydrogen negative ion H~ , was sug- 
gested in the fifties ^^ and is now well established exper- 
imentally. Negatively charged donors in bulk Si were 
first detected by photoconductivity measurements J^ The 
binding energies of D~ donors, defined as the energy re- 
quired to remove one electron from the ion (D^ — > D^+ 
free-electron) Eg = Edo — Ejj- , are found experimen- 
tally to be small {E^' ~ 1.7 meV for P and - 2.05 
meV for As) compared to the binding energies of the 
first electron Eg (45 meV for P and 54 meV for As). 
For zero applied magnetic fields, no excited bound states 
of I?^ in bulk semiconductorsi^ or superlatticesi^ are 
found, similar to H~ which has only one bound state in 
three-dimensions as shown in Refs. 20.21. 

A relevant characteristic of negatively charged donors 
is their charging energy, U — Ejj- — 2Euo, which gives 
the energy required to add a second electron to a neutral 
donor. This extra energy is due to the Coulomb repulsion 
between the two bound electrons, and does not contribute 
in one electron systems, as D^. The measured values in 



bulk Si are U 
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meV for P. 

From the stability diagrams obtained from transport 
spectroscopy measurements we observe that the charging 
energy of As dopants in nanoscale Si devices (FinFETs) is 
strongly reduced compared to the well known bulk value. 
By using a variational approach within the single- valley 
effective mass approximation, we find that this decrease 
of the charging energy may be attributed to modifica- 
tions on the bare insulator screening by the presence of a 
nearby metallic layer. For the same reason, we also find 
theoretically that it may be possible to have a. D^ bound 
excited state. 

This paper is organized as follows. In Sec. [TTl we in- 
troduce the formalism for a donor in the bulk in analogy 
with the hydrogen atom problem. In Sec. IIIIl we study 



the problem of a donor close to an interface within a flat 
band condition. We show experimental results for the 
charging energy and compare them with our theoretical 
estimations. We also calculate the binding energy of a 
D^ triplet first excited state. In Sec. llVl we present dis- 
cussions including: (i) assessment of the limitations in 
our theoretical approach, (ii) considerations about the 
modifications of the screening in nanoscale devices, (iii) 
the implications of our results in quantum device appli- 
cations, and, finally, we also summarize our main conclu- 
sions. 



II. DONORS IN BULK SILICON 

A simple estimate for the binding energies of both D^ 
and D~ in bulk Si can be obtained using the analogy 
between the hydrogen atom H and shallow donor states 
in semiconductors. The Hamiltonian for one electron in 
the field of a nucleus with charge -j-e and infinite mass 



is, in effective units of length as 



Ry 



h /rriee and energy 
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with T{r) = — V^. The ground state is 
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with Bohr radius a — las and energy Eh = —IRy. 
This corresponds to one electron in the Is orbital. 

For negatively charged hydrogen (H^) the two elec- 
trons Hamiltonian is 



^Buik = h{ri) + h{r2) H , 

ri2 



(3) 



where the last term gives the electron-electron interac- 
tion (ri2 = |ri — r2|)- As an approximation to the ground 
state, we use a relatively simple variational two particles 
wave-function for the spatial part, a symmetrized combi- 
nation of Is atomic orbitals as given in Eq. [2l since the 
spin part is a singlet, 

|ls,ls,s) = [(/)(ri,a)0(r2,6)-f (/)(ri,6)(/)(r2,a)] . (4) 

The resulting energy is E^ = —l.Q27Ry with a = 
0.963 Ob and h — 3.534 a^ (binding energy E^ = 
O.Q27Ry).^^ Here we may interpret a as the radius of 
the inner orbital and b of the outer orbital. This approx- 
imation for the wave-function correctly gives a bound 
state for H~ but it underestimates the binding energy 
with respect to the value E§ = 0.0555 Ry, obtained 
with variational wave-functions with a larger number of 
parameters, thus closer to the 'exact' value. ^^ 

Assuming an isotropic single-valley conduction band 
in bulk Si the calculation of the D" and D~ energies 
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TABLE I: Bulk values of energies and orbital radii for the 
ground state of neutral and negatively charged donors within 
our approximation (see text for discussion). Effective units 
for Si are a* — 2.14 nm and Ry* = 31.2 meV. 
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FIG. 1: (Color online) Schematic representation of a nega- 
tively charged donor in Si (solid circles) located a distance 
d from an interface. The open circles in the barrier (left) 
represent the image charges. The sign and magnitude of 
these charges depend on the relation between the dielectric 
constants of Si and the barrier given by Q = (ebarrior — 
esi)/(ebarrier + ESi)- For the electrous, Q < corresponds to 
repulsive electron image potentials and a positive donor im- 
age potential (opposite signs of potentials and image charges 
for Q>0, seeEq.El). 



reduces to the case of H just described. Within this ap- 
proximation, an estimation for Eg can then be obtained 
by considering an effective rydberg Ry* — m*e^/2egfi^ 
with an isotropic effective mass (we use esi = 11.4). We 
choose m* = 0.29819 We so that the ground state energy 
for a neutral donor is the same as given by an anisotropic 
wave-function in bulk: within a single valley approxima- 
tion Eg = —IRy* = —31.2 meV and its effective Bohr 
radius is o = la* with a* = Ti esi/m*e^ = 2.14 nm. In 
this approximation. Eg — 0.84 meV. In the same way, 
an estimation for the charging energy can be made for 
donors in Si: U = 0.973Ry* = 30.35 meV.^^ 

Even though the trial wave function in Eq. ^ underes- 
timates the binding energy, we adopt it here for simplic- 
ity, in particular to allow performing in a reasonably sim- 
ple way the calculations for a negatively charged donor 
close to an interface reported below. In the same way, we 
do not introduce the multivalley structure of the conduc- 
tion band of Si. The approximations proposed here lead 
to qualitative estimates and establish general trends for 
the effects of an interface on a donor energy spectrum. 
The limitations and consequences of our approach are 
discussed in Sec. IIVI 
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FIG. 2: (Color online) Energy of the neutral donor versus 
its distance d from an interface for different values of Q = 

(.^barrier ^Si j/ l^^barricr i ^Si}' 



III. DONORS CLOSE TO AN INTERFACE 
A. D° and D^ ground states 

We consider now a donor (at z ~ 0) close to an inter- 
face (at z — —d) (see Fig.[T]). Assuming that the interface 
produces an infinite barrier potential, we adopt varia- 
tional wave-functions with the same form as in Eqs. ([2]) 
and ^ multiplied by linear factors {zi + d) {i = 1,2) 
which guarantee that each orbital goes to zero at the in- 
terface. We further characterize the Si interface with a 
different material by including charge image terms in the 
Hamiltonian. 

Before discussing the ionized donor D^ , we briefly 
present results for the neutral donor D'^ which are in- 
volved in defining donor binding and charging energies. 
For this case, the Hamiltonian is 



H{ri) = h{ri) + himngcsin) 
with h{ri) as in Eq. [T]and 
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h; 



images 



(^i) = - 



Q 



2Q 



2(^1 + d) ^/xf+W+W+W ' 



(6) 

where Q = (ebanicr ^ eSi)/(ebarrior + <^Si)- ebarrior is the 

dielectric constant of the barrier material. The first term 
in ^images is the interaction of the electron with its own 
image, and the second is the interaction of the electron 
with the donor's image. If the barrier is a thick insulator, 
for example Si02 with dielectric constant esi02 = 3.8, 
Q < {Q = —0.5 in this case). In actual devices, the 
barrier is composed of a thin insulator (usually Si02), 
which prevents charge leakage, plus metallic electrodes 
which control transport and charge in the semiconductor. 
This composite heterostructure may effectively behave as 
a barrier with an effective dielectric constant larger than 
Si, since e,nctai — > oo, leading to an effective Q > 0. 
Depending on the sign of Q, the net image potentials 
will be repulsive or attractive, which may strongly affect 
the binding energies of donors at a short distance d from 
the interface. 

Using a trial wave- function (pDO ex e~^'°'{z + d), most 
of the integrals involved in the variational calculation of 



E'^ can be performed analytically. E^ is shown in 
Fig. [2] for different values of Q and compare very well with 
the energy calculated by MacMillen and Landman^ for 
Q — —1 with a much more complex trial wave- function. 
The main effect of the interface is to reduce the binding 
energy when the donor is located at very small distances 
d. For Q < (corresponding to insulating barriers with 
a dielectric constant smaller than that of Si) , the energy 
has a shallow minimum for d ~ 8a*. This minimum arises 
because the donor image attractive potential enhances 
the binding energy but, as d gets smaller, the fact that 
the electron's wave- function is constrained to z > ~d 
dominates, leading to a strong decrease in the binding 
energy,^ Q — corresponds to ignoring the images. Q = 
1 would correspond to having a metal at the interface 
with an infinitesimal insulating barrier at the interface to 
prevent leackage of the wave-function into the metal.— 
We show results for Q = 0.5 as an effective value to take 
account of a realistic barrier composed of a thin (but 
finite) insulator plus a metal. The bulk limit E — —1 Ry* 
is reached at long distances for all values of Q. 

Adding a second electron to a donor requires the inclu- 
sion of the electron-electron interaction terms. The neg- 
ative donor Hamiltonian parameters are schematically 
presented in Fig. [TJ and the total two electrons Hamilto- 
nian is 



H = Hin) + H{r2) 



ri2 



4Q 



,(7) 



y/{xi - X2f + (yi - 2/2)2 _^ (2d + zi + za)- 

where H{7\) includes the one-particle images (Eq.[5]) and 
the last term is the interaction between each electron and 
the other electron's image. 

In Figs. |3] and |4l we plot E^ and the binding en- 
ergy Eg — E^ — E^ assuming a trial wave-function 
ex [(f>{ri,a)(f>{r2,b) + (j){ri,b)(f>{r2,a)] (zi +d){z2 + d) with 
variational parameters a and b, for Q ~ —0.5 and Q — 0.5 
respectively. The radius of the inner orbital is a ^ la* 
while 6, the radius of the outer orbital, depends very 
strongly on Q and d and is shown in Figs.|3Uc) andUl^c). 
We have done calculations for several values of Q, rang- 
ing from Q = +\ to Q = —\. The general trends and 
qualitative behavior of the calculated quantities versus 
distance d are the same for all Q > (effective barrier 
dominated by the metallic character of the interface ma- 
terials), which differ from the also general behavior of 
Q < (effective barrier dominated by the insulator ma- 
terial). For Q < (illustrated for the particular case of 
Q — —0.5 in Fig. [3)), D^ is not bound for small d (for 
d < 4a* in the case oi Q — —0.5). For larger d's, the 
binding energy is slightly enhanced from the bulk value. 
The radius of the outer orbital b is very close to the bulk 
value for d > 4 a*. For Q > (illustrated by Q = 0.5 
in Fig. 2]), D^ is bound at all distances d, though the 
binding energy is smaller than in bulk. The radius of the 
outer orbital b is very large and increases linearly with 
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FIG. 3; (Color online) Results for Q = —0.5. (a) Energy for 
a neutral donor D", and for the ground D~ |ls, Is, s) and first 
excited Z)~|ls,2s,f) negatively charged donor, (b) Binding 
energies of the D~ states, (c) Value of the variational param- 
eter h for the D~ ground state. For d < 4, D~ |ls, Is, s) is not 
stable and the energy is minimized with 6 -^ oo. Bulk values 
are represented by short line segments on the right. 



d up to ^crossover ~ 14.5a* [see Fig. S] (c)]. For larger 
d, h is suddenly reduced to its bulk value. This abrupt 
behavior of the h that minimizes the energy is due to two 
local minima in the energy versus h: for d < dcrossovor 
the absolute minimum corresponds to a very large (but 
finite) orbital radius b while for d > d. 



'crossover 



the abso- 
lute minimum crosses over to the other local minimum, 
at 6 ~^ fcbuik- As d increases from the smallest values and 
b increases up to the discontinuous drop, a "kink" in the 
D~ binding energy is obtained at the crossover point [see 
Fig. U (b)], changing its behavior from a decreasing to 
an increasing dependence on d towards the bulk value as 
d — >■ oo. 



B. Charging energy: experimental results 

The charging energy of shallow dopants can be ob- 
tained by using the combined results of photoconductiv- 
ity experiments to determine the D~ binding energy-^ 
and direct optical spectroscopy to determine the binding 
of the -D'^ state.2S It was shown recently that the charging 
energy in nanostructures can be obtained directly from 
charge transport spectroscopy at low temperature.^^ Sin- 
gle dopants can be accessed electronically at low temper- 
ature in deterministically doped silicon/silicon-dioxide 
heterostructures^^ and in small silicon nanowire field ef- 
fect transistors (FinFETs), where the dopants are posi- 
tioned randomly in the channel J^ii^i^ Here we will fo- 
cus in particular on data obtained using the latter struc- 
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FIG. 4: (Color online) Same as Fig.[3]for Q = 0.5. 



tures.— i^l 

FinFET devices in which single dopant transport have 
been observed typically consist of crystalline silicon wire 
channels with large patterned contacts fabricated on 
silicon-on-insulator. Details of the fabrication can be 
found in Ref. [l^. In this kind of samples, few dopants 
may diffuse from the source/drain contacts into the chan- 
nel during the fabrication modifing the device character- 
istics both at rooroi^ and low temperaturesii^ii^'21 In 
some cases, subthreshold transport is dominated by a 
single dopant,^ 

Low temperature transport spectroscopy relies on the 
presence of efficient Coulomb blockade with approxi- 
mately zero current in the blocked region. This requires 
the thermal energy of the electrons, fc^T, to be much 
smaller than C/, a requirement that is typically satisfied 
for shallow dopants in silicon at liquid helium temper- 
ature and below, i.e. < A.2K. At these temperatures 
the current is blocked in a diamond-shaped region in a 
stability diagram, a color-scale plot of the current - or 
differential conductance dl/dVb - as a function of the 
source/drain, 14, and gate voltage, Vg. 

In Fig. [51 the stability diagram of a FinFET with only 
one As dopant in the conduction channel is shown. At 
small bias voltage (eVb <^ ksT), increasing the voltage 
on the gate effectively lowers the potential of the donor 
such that the different donor charge states can become 
degenerate with respect to the chemical potentials in the 
source and drain contacts and current can flow. The dif- 
ference in gate voltage between the D^/D" and D"/D~ 
degeneracy points (related to the charging energy) de- 
pends, usually in good approximation, linearly on the 
gate voltage times a constant capacitive coupling to the 
donorJ^ Generally a more accurate and direct way to 
determine the charging energy is to determine the bias 
voltage at which the Coulomb blockade for a given charge 



Differential Conductance [|iS] 



> 

E 




220 240 260 280 300 320 340 360 380 400 420 
Vg[mV] 

FIG. 5: (Color online) DifTerential conductance stability di- 
agram showing the transport characteristics of a single As 
donor in a FinFET device. ~ The differential conductance is 
obtained by a numerical differentiation of the current with re- 
spect to Vb at a temperature of 0.3 K. Extracting the charging 
energy from the stability diagram can be done by determin- 
ing the gate voltage for which Coulomb blockade of a given 
charge state (the D" charge state in this case) is lifted for all 
Vg. The transition point is indicated by the horizontal arrow, 
leading to a charging energy f/ = 36 meV, as given by the 
vertical double-arrow. The inset shows the electrical circuit 
used for the measurement. 



state is lifted for all gate voltages, indicated by the hor- 
izontal arrow in Fig. [S] This method is especially useful 
when there is efficient Coulomb blockade. For the partic- 
ular sample shown in Fig. [SI f/ = 36 meV. This is similar 
to other reported values in the literatur©iiii^i21 ranging 
from ^ 26 to ^ 36 meV. There is therefore a strong 
reduction in the charging energy compared to the bulk 
value t/buik — 52 mcV. The ratio between the observed 
and the bulk value is ^ 0.6 — 0.7. 

Theoretically, we can extract the charging energy from 
the results in Figs. |3] and S) The results are shown as a 
function of d for Q — —0.5, 0, and 0.5 in Fig. [B] A 
reduction of the charging energy U of the order of the 
one observed occurs at d ^ 2a* for 0.1 < Q < 1 (only 
Q — 0.5 is shown in the figure). Therefore, the exper- 
imentally observed behavior of U is consistent with a 
predominant influence of the metallic gates material in 
the D~ energetics. On the other hand, for Q < 0, U 
is slightly enhanced as d decreases and, for the smallest 
values of d considered, the outer orbital is not bound. At 
very short distances d, the difference in behavior between 
the insulating barrier {Q < 0) and the barrier with more 
metallic character (Q > 0) is in the interaction between 
each electron and the other electron's image, which is re- 
pulsive in the former case and attractive in the latter. 
Although this interaction is small, it is critical to lead to 
a bound D~ for Q > and an unbound D^ for Q < 
at very short d. 
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FIG. 6: (Color online) Charging energy U of the D~ ground 
state for three different values of Q. For Q < 0, the charging 
energy is nearly constant with d. For these cases, the nega- 
tively charged donor is not bound for small d. For Q > the 
charging energy decreases as the donor gets closer to the inter- 
face, at relatively small distances d. The latter is consistent 
with the experimental observation. 



D first excited state. 



It is well established that in 3 dimensions (with no 
magnetic field applied) there is only one bound state 
of D~ J^^ Motivated by the significant changes in the 
ground state energy produced by nearby interfaces, we 
explore the possibility of having a bound excited state in 
a double-charged single donor. Like helium, we expect 
the D^ first excited state to consist of promoting one Is 
electron to the 2s orbital. The spin triplet |ls, 2s, t) state 
(which is orthogonal to the singlet ground state) has a 
lower energy than |ls,2s,s)i^ As a trial wave-function 
for I Is, 2s, t) we use the antisymmetrized product of the 
two orbitals Is and 2s and multiply by (zi -I- d){z2 + d) 
to fulfill the boundary condition, namely. 



^1 



s,2s,t 



= N 



\2b J 



- i -e 



\2b J\ 



X {zi+d){z2+d) 



(8) 



with a and b variational parameters and A'^ a normaliza- 
tion factor. Note that, for a particular value of b, the 
outer electron in a 2s orbital would have a larger effec- 
tive orbital radius than in a Is orbital due to the different 
form of the radial part. 

For Q < 0, the outer orbital is not bound and the 
energy reduces to that of D^ (see Fig. [3l). Surprisingly, 
for Q > the |ls, 2s, t) state is bound and, as d increases, 
tends very slowly to the D'^ energy as shown in Fig. H) 
Moreover, its binding energy is roughly the same as the 
ground state |ls. Is, s) for d < 15a*, another unexpected 
result. The existence of a bound D^ triplet state opens 
the possibility of performing coherent rotations involving 
this state and the nearby singlet ground state. 



IV. DISCUSSIONS AND CONCLUSIONS 

Our model for D^ centers involves a number of sim- 
plifications: (i) the mass anisotropy is not included; (ii) 
the multivalley structure of the conduction band of Si is 
not considered; (iii) correlation terms in the trial wave- 
function are neglected. These assumptions aim to de- 
crease the number of variational parameters while allow- 
ing many of the integrals to be solved analytically. 

Qualitatively, regarding assumption (i), it has been 
shown that the mass anisotropy inclusion gives an in- 
crease of the binding energy for both D° and D^ (see 
Ref. \2§) ; regarding (ii) , inclusion of the multivalley struc- 
ture of the conduction together with the anisotropy of the 
mass would lead to an enhancement of the binding energy 
of D~ due to the possibility of having intervalley config- 
urations in which the electrons occupy valleys in 'per- 
pendicular' orientations, (with perpendicularly oblated 
wave- functions), thus leading to a strong reduction of 
the electron-electron repulsive interaction<22i2£ Regard- 
ing point (iii), more general trial wave- functions for D~ 
have been proposed in the literature. For example, the 
one suggested by Chandrasekar models correlation effects 
by multiplying Eq. |3]by a factor, (1 -I- Cri2), ~ where C 
is an additional variational parameter. In the bulk, the 
effect of this correlation factor is to increase the bind- 
ing energy of D~ from 0.027 Ry* if C = (our case) to 
0.0518 i?j/*<2^ We conclude that all three simplifications 
assumed in our model lead to an underestimation of the 
binding energy of D~ , thus, the values reported here are 
to be taken as lower bounds for it. 

As compared to experiments, an important difference 
with respect to the theory is that we are assuming a flat- 
band condition while the actual devices have a built-in 
electric field due to band-bending at the interface be- 
tween the gate oxide and the p-doped channel>21yS If 
an electric field were included, the electron would feel a 
stronger binding potential (which results from the addi- 
tion of the donor potential and the triangular potential 
well formed at the interface) leading to an enhancement 
of the binding energy of Z?" and D^ (with an expected 
strong decrease of the electron-electron interaction in this 
case for configurations with one electron bound to the 
donor at z = and the other pulled to the interface at 
z = -d). 

The presented results are dominated by the presence 
of a barrier, which constrains the electron to the z > —d 
region, and the modification of the screening due to the 
charge induced at the interface, a consequence of the di- 
electric mismatch between Si and the barrier material. 
This is included by means of image charges. Effects of 
quantum confinement and dielectric confinement^ii^ are 
not considered here: we believe these are not relevant 
in the FinFETs under study. Although the conduction 
channel is very narrow (4 nm^ )^ the full cross section of 
the Si wire is various tens of nm and quantum and di- 
electric confinement is expected to be effective for typical 
device sizes under 10 nm. Both quantum and dielectric 



confinement lead to an enhancement of the binding en- 
ergy with respect to the bulk, which is the opposite to 
what we obtain for small d. 

Neutral double donors in Si, such as Te or Se, have been 
proposed for spin readout via spin-to-charge conversion^i 
and for spin coherence time measurements.^^ The nega- 
tive donor D~ also constitutes a two-electron system, 
shallower than Te and Se. In this context, investigation 
of the properties oi D~ shallow donors in Si affecting 
quantum operations as, for example, their adequacy for 
implementing spin measurement via spin-to-charge con- 
version mechanism, ^4'-?§ deserve special attention. Our 
theoretical study indicates that, very near an interface 
(for d < 4a*), the stability of D~ against dissociation 
requires architectures that yield effective dielectric mis- 
match Q > 0, a requirement for any device involving 
operations or gates based on D~ bound states. 

In conclusion, we have presented a comprehensive 
study of the effects of interface dielectric mismatch in 
the charging energy of nearby negatively charged donors 
in Si. In our study, the theoretical treatment is based on 
a single- valley effective mass formalism, while transport 
spectroscopy experiments were carried out in FinFET de- 
vices. The experiments reveal a strong reduction on the 
charging energy of isolated As dopants in FinFETs as 
compared to the bulk values. Calculations present, be- 
sides the charging energy, the binding energy of donor in 
three different charge states as a function of the distance 
between the donor and an interface with a barrier. The 
boundary problem is solved by including the charge im- 
ages whose signs depend on the difference between the 
dielectric constant of Si and that of the barrier material 
[the dielectric mismatch, quantified by the parameter Q 
defined below Eq. (|6])]. 

Typically, thin insulating layers separate the Si chan- 
nel, where the dopants are located, from metallic gates 
needed to control the electric fields applied to the de- 
vice. This heterostructured barrier leads to an effective 
screening with predominance of the metallic components, 
if compared to a purely Si02 thick layer, for which Q < 0. 
Assuming a barrier material with an effective dielectric 
constant larger than that of Si (in particular, Q = 0.5 
corresponds to ebarrier = 3esi), we obtain a reduction of 
the charging energy U relative to Uhuik at small d, consis- 
tent with the experimental observation. We did not at- 
tempt quantitative agreement between presented values 
here, but merely to reproduce the right trends and clar- 
ify the underling physics. It is clear from our results that 
more elaborate theoretical work on interface effects in 
donors, beyond the simplifying assumptions here, should 
take into account the effective screening parameter as a 
combined effect of the nearby barrier material and the 
adjacent metallic electrodes. From our calculations and 
experimental results, we conclude that the presence of 
metallic gates tend to increase fbardo™ above esi, leading 
to Q > and reducing the charging energies. 
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